Method for programming repeating motion of redundant robotic arm

ABSTRACT

A method is presented for programming a repeating motion of a redundant robotic arm on the basis of a variable parameter convergence differential neural network. The method may include establishing an inverse kinematics equation, creating an inverse kinematics problem, introducing a repeating motion indicator, converting a time-varying convex quadratic programming problem into a time-varying matrix equation, and integrating an optimal solution to obtain an optimal solution of a joint angle. The use of the variable parameter convergence differential neural network to solve the repeating redundant mechanical motion has the advantages of high computational efficiency, high real-time performance, and enhanced robot arm robustness.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority to China Patent Application No. CN 201710188430.4 filed Mar. 27, 2017, and International Patent Application No. PCT/CN2017/111107 filed Nov. 15, 2017, both of which are hereby incorporated by reference in their entirety.

FIELD OF THE DISCLOSURE

The present disclosure relates to the field of redundant robotic arms, and in particular to a method for programming and operating a repeating motion of a redundant robotic arm on the basis of a variable parameter convergence differential neural network.

BACKGROUND OF THE DISCLOSURE

The term “redundant robotic arm” means that the number of degrees of freedom of the robotic arm is greater than the number of degrees of freedom required to complete a task, so that when completing the main tasks by an end effector, the redundant robotic arm can also perform additional tasks such as avoiding an obstacle avoidance, a shutdown extreme position, and a singular state of the robotic arm due to the more degrees of freedom. The term “repeating motion” means that after the end of the robotic arm completes a cycle of actions, all the joints thereof can return to the initial positions, not just the end of the robotic arm returns back to the initial position. In the automated industrial production, the robotic arm usually needs to carry out mass production activities, if the robotic arm completes a non-repeating motion, that is, the initial state of each cycle of motions is different, an error will occur, and after the error has been accumulated to a certain extent, additional reset operation of the robotic arm is required, and the production efficiency will be greatly reduced. Therefore, it is meaningful to study the repeating motion of the redundant robotic arm.

The traditional solution to the inverse kinematics problem of the redundant robotic arm is based on a pseudo-inverse method, this method is computationally intensive, has poor real-time performance, and considers a single problem constraint, which is greatly restricted in the actual application of the robotic arm. A quadratic programming scheme has been proposed to solve the repeating motion of a redundant robotic arm, which is divided into a numerical method solver and a neural network solver. Compared with the numerical method solver, the neural network solver has the advantages of being more efficient and having better real-time performance.

OVERVIEW OF THE DISCLOSURE

In view of the deficiencies of the prior methodologies, an object of the present disclosure is to provide a method for programming a repeating motion of a redundant robotic arm on the basis of a variable parameter convergence differential neural network, which has a higher calculation accuracy and a better robustness compared with the classical recurrent neural network solver to solve the repeating motion problem of a redundant robotic arm.

While other objects of this disclosure are contemplated, for example purposes only, an object of the present disclosure may be achieved by means of the following exemplary technical solution:

A method for programming a repeating motion of a redundant robotic arm on the basis of a variable parameter convergence differential neural network, the method including the steps of:

1) establishing an inverse kinematics equation of the redundant robotic arm at a velocity level by means of a track of an end of the redundant robotic arm;

2) designing an inverse kinematics problem in step 1 as a time-varying convex quadratic programming problem constrained by an equality;

3) introducing a repeating motion indicator into the time-varying convex quadratic programming problem of step 2;

4) converting the time-varying convex quadratic programming problem, into which the repeating motion indicator is introduced in step 3, into a time-varying matrix equation by using a Lagrangian function;

5) solving the time-varying matrix equation of step 4 by means of the variable parameter convergence differential neural network; and

6) integrating an optimal solution, obtained in step 5, of the redundant robotic arm at the velocity level to obtain an optimal solution of a joint angle.

Further, in step 1, the inverse kinematics equation of the redundant robotic arm may be expressed as: f(θ)=r

where r is a desired track of the end of the redundant robotic arm, and f(•) is a nonlinear equation of the joint angle of the redundant robotic arm, and the inverse kinematics equation of the redundant robotic arm at the velocity level is obtained by deriving the two sides of the equation with respect to time: J(θ){dot over (θ)}={dot over (r)}

where J(θ)∈R^(m×n) is an m×n-dimensional matrix on the real number field, J(θ) is a Jacobian matrix of the redundant robotic arm, n is the number of the degrees of freedom of the robotic arm, and m is the number of the spatial dimensions of the track of the end of the robotic arm, and {dot over (θ)} and {dot over (r)} are respectively the derivatives of the joint angle and the track of the end of the redundant robotic arm with respect to time.

Further, in step 2, the specific formula for designing the inverse kinematics problem in step 1 as a time-varying convex quadratic programming problem constrained by an equality may be:

$\min\left( {\frac{x^{T}{Wx}}{2} + {c^{T}x}} \right)$ s.t.  (J(θ)x = b)

where x={dot over (θ)}, b={dot over (r)}. W=I represents an identity matrix, J(θ) is the Jacobian matrix of the redundant robotic arm, and c is a performance indicator coefficient.

Further, the repeating motion indicator in step 3 is obtainable by means of the performance indicator coefficient c, which is designed as c=ζ(θ(t)−θ(0)), where ζ represents the response coefficient to a joint offset, and θ(t) and θ(0) respectively represent a joint state during the movement of the robotic arm and an initial joint state.

Further, in step 4, the Lagrangian function may be constructed as:

${L\left( {x,\lambda,t} \right)} = {\frac{x^{T}{Wx}}{2} + {c^{T}x} + {\lambda\left( {{{J(\theta)}x} - b} \right)}}$

where λ is a Lagrangian multiplier, and the partial derivatives of the Lagrangian function respectively with respect to x and λ are obtained:

$\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial x} = {{{Wx} + c + {{J^{T}(\theta)}\lambda}} = 0}$ $\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial\lambda} = {{{{J(\theta)}x} - b} = 0}$

the above system of equations can be expressed as the following time-varying matrix equation:

     Qy = u ${{{where}\mspace{14mu} Q} = {\begin{bmatrix} W & {J^{T}(\theta)} \\ {J(\theta)} & 0 \end{bmatrix} \in E^{{({n + m})} \times {({n + m})}}}},{y = {\begin{bmatrix} x \\ \lambda \end{bmatrix} \in R^{({n + m})}}},{u = {\begin{bmatrix} {- c} \\ b \end{bmatrix} \in R^{({n + m})}}}$

Further, the specific process of step 5 may be that the error of the time-varying matrix equation converges to zero based on the variable parameter convergence differential neural network, firstly, an error function is constructed as: ε(t)=Qy−u

where ε(t) represents the error of the time-varying matrix equation, and then based on a neurodynamic method, the error may be designed to converge to zero in the following way, the specific formula is:

$\frac{d\;{ɛ(t)}}{dt} = {{- \gamma}\;{\exp(t)}{\Phi\left( {ɛ(t)} \right)}}$

where γ may be a parameter for adjusting a convergence rate, Φ(•) is an activation function, and the error function is substituted into the above formula to obtain a variable parameter convergence differential neural network solver, namely: Q{dot over (y)}=−{dot over (Q)}y−γ exp(t)Φ(Qy−u)+{dot over (u)}

In this way, the variable parameter convergence differential neural network solver obtains an optimal solution y* of the time-varying matrix equation, and the first n terms thereof are the optimal solution x* of the time-varying convex quadratic programming problem in step 2 i.e., the optimal solution of a joint angular velocity.

Further, the optimal solution θ* of the joint angle in step 6 is obtained by integrating the optimal solution of the time-varying convex quadratic programming problem, i.e. the optimal solution x* of the joint angular velocity.

As compared with prior systems, the present disclosure has the following advantages and beneficial effects:

The present disclosure may solve the problem of non-repeating motion of the redundant robotic arm by means of the quadratic programming scheme, which can consider various constraint conditions, and has a good real-time performance as compared with the traditional pseudo-inverse based method.

The present disclosure may use a neural network solver, and the calculation speed is faster and the efficiency is higher as compared with the numerical method solver.

The present disclosure may use a novel variable parameter convergence differential neural network solver, which has a faster convergence speed, a higher calculation precision and a better robustness as compared with the classical recurrent neural network solver.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flowchart of a method for programming a repeating motion of a redundant robotic arm of an embodiment of the present disclosure;

FIG. 2(a) is a schematic diagram showing the redundant robotic arm of the embodiment of the present disclosure performing a non-repeating motion; and

FIG. 2(b) is a schematic diagram showing the redundant robotic arm of the embodiment of the present disclosure performing a repeating motion.

DETAILED DESCRIPTION OF THE EMBODIMENTS

The present disclosure will be further described in detail below in connection with embodiments and the accompanying drawings, but embodiments of the present disclosure are not limited thereto.

An exemplary embodiment provides a method for programming a repeating motion of a redundant robotic arm 10 on the basis of a variable parameter convergence differential neural network, the flow chart thereof is shown in FIG. 1, and is composed of three parts: presentation of a problem 12, conversion of the problem 14 and creating a solution to the problem 16. Firstly, the inverse kinematics equation at a velocity level is established according to the desired track of an end of the redundant robotic arm and the Jacobian matrix. The repeating motion of the redundant robotic arm is designed as a time-varying convex quadratic programming problem constrained by an equality. The quadratic programming problem is converted into a matrix equation problem by means of a Lagrangian function. And finally, the matrix equation is solved by means of the variable parameter convergence differential neural network. The following steps are contemplated. It will be appreciated that more or fewer steps could be deployed:

1) establishing an inverse kinematics equation 20 of the redundant robotic arm 10 at a velocity level by means of a track 18 of an end of the redundant robotic arm. In this step, the inverse kinematics equation 20 of the redundant robotic arm is expressed as: f(θ)=r

where r is a desired track of the end of the redundant robotic arm, and f(•) is a nonlinear equation of the joint angle of the redundant robotic arm, which is determined by the structure of the robotic arm 10, a Kinova Jaco six-axis robotic arm is simulated in this embodiment, and the inverse kinematics equation 20 of the redundant robotic arm at the velocity level is obtained by deriving the two sides of the equation with respect to time: J(θ){dot over (θ)}={dot over (r)}

Here, J(θ)∈R^(m×n) is an m×n-dimensional matrix on the real number field, J(θ) is a Jacobian matrix of the redundant robotic arm, n is the number of the degrees of freedom of the robotic arm, and m is the number of the spatial dimensions of the track of the end of the robotic arm, and {dot over (θ)} and {dot over (r)} are respectively the derivatives of the joint angle and the track of the end of the redundant robotic arm 10 with respect to time;

2) Designing or creating an inverse kinematics problem in step 1 as a time-varying convex quadratic programming problem constrained by an equality.

An exemplary formula is:

$\min\left( {\frac{x^{T}{Wx}}{2} + {c^{T}x}} \right)$ s.t.  (J(θ)x = b)

where x={dot over (θ)}, b={dot over (r)}. W=I represents an identity matrix, J(θ) is the Jacobian matrix 24 of the redundant robotic arm 10, and c is a performance indicator coefficient.

3) Introducing a repeating motion indicator 26 into the time-varying convex quadratic programming problem of step 2. The repeating motion indicator in this step is obtainable by means of the performance indicator coefficient c, which is designed as c=ζ(θ(t)−θ(0)), where ζ represents the response coefficient to a joint offset, and in this embodiment ζ=5; and θ(t) and θ(0) respectively represent a joint state during c the movement of the robotic arm and an initial joint state. When the repeating motion indicator c is not considered, ζ=0, the simulation result is as shown in FIG. 2(a). In this case, after the robotic arm 10 completes a cycle of motions, it is impossible to ensure that the end position 11 and the initial position 13 of each joint of the robotic arm 10 are consistent; and when the repeating motion indicator is considered, ζ=5, the simulation result is as shown in FIG. 2(b), the end position 11 of each joint of the robotic arm 10 is consistent with the initial position 13 thereof. This arrangement is preferred.

4) Converting the time-varying convex quadratic programming problem, into which the repeating motion indicator is introduced in step 3, into a time-varying matrix equation 30 by using a Lagrangian function 28.

An example of a process of this step is that the Lagrangian function 28 is constructed as:

${L\left( {x,\lambda,t} \right)} = {\frac{x^{T}{Wx}}{2} + {c^{T}x} + {\lambda\left( {{{J(\theta)}x} - b} \right)}}$

where λ is a Lagrangian multiplier, and the partial derivatives of the Lagrangian function 28 respectively with respect to x and λ are obtained:

$\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial x} = {{{Wx} + c + {{J^{T}(\theta)}\lambda}} = 0}$ $\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial\lambda} = {{{{J(\theta)}x} - b} = 0}$

The above system of equations can be expressed as the following time-varying matrix equation 30:

     Qy = u ${{{where}\mspace{14mu} Q} = {\begin{bmatrix} W & {J^{T}(\theta)} \\ {J(\theta)} & 0 \end{bmatrix} \in E^{{({n + m})} \times {({n + m})}}}},{y = {\begin{bmatrix} x \\ \lambda \end{bmatrix} \in R^{({n + m})}}},{u = {\begin{bmatrix} {- c} \\ b \end{bmatrix} \in {R^{({n + m})}.}}}$

5) Solving the time-varying matrix equation 30 of step 4 by means of the variable parameter convergence differential neural network 32; and the specific process of this step is that the error of the time-varying matrix equation converges to zero based on the variable parameter convergence differential neural network. Firstly, an error function is constructed as: ε(t)=Qy−u

where ε(t) represents the error of the time-varying matrix equation, and then based on a neurodynamic method, the error is designed to converge to zero in the following way, the specific formula is:

$\frac{d\;{ɛ(t)}}{dt} = {{- \gamma}\;{\exp(t)}{\Phi\left( {ɛ(t)} \right)}}$

where γ is a parameter for adjusting a convergence rate, Φ(•) is an activation function, and the error function is substituted into the above formula to obtain a variable parameter convergence differential neural network solver 32, namely: Q{dot over (y)}=−{dot over (Q)}y−γ exp(t)Φ(Qy−u)+{dot over (u)}

In this way, the variable parameter convergence differential neural network solver 32 obtains an optimal solution y* of the time-varying matrix equation, and the first n terms thereof are the optimal solution x* of the time-varying convex quadratic programming problem in step 2 i.e., the optimal solution of a joint angular velocity.

6) Integrating an optimal solution, obtained in step 5, of the redundant robotic arm 10 at the velocity level to obtain an optimal solution of a joint angle 34.

The optimal solution θ* of the joint angle 34 in this step is obtained by integrating the optimal solution of the time-varying convex quadratic programming problem, i.e. the optimal solution x* of the joint angular velocity.

The specific process of this step is that by means of the variable parameter convergence differential neural network solver in step 5, the optimal solution y* can be obtained, and the first n terms thereof are the optimal solution x* of the time-varying convex quadratic programming problem in step 2 i.e., the optimal solution of a joint angular velocity, which can be integrated to obtain the optimal solution θ* of the joint angle of the redundant robotic arm.

The foregoing description is merely illustrative of the disclosed embodiments of the present disclosure, but the scope of protection of the present disclosure is not limited thereto. Equivalent replacements or modifications made to the technical solutions and the inventive concept of the present disclosure by a person skilled in the art within the scope of the disclosure of the present disclosure fall into the scope of protection of the present disclosure. 

The invention claimed is:
 1. A method for performing a repeating motion of a robotic arm on the basis of a variable parameter convergence differential neural network, comprising: providing a robotic arm, the robotic arm having an end and a plurality of joints; providing a system for operating the robotic arm by the following steps: a) establishing an inverse kinematics equation of the robotic arm at a velocity level by tracking the end of the robotic arm; b) generating an inverse kinematics problem as a time-varying convex quadratic programming problem constrained by an equality; c) obtaining an initial joint state θ(0) of the robotic arm and a joint state θ(t) during the motion of the robotic arm, and introducing a repeating motion indicator into the time-varying convex quadratic programming problem; wherein the repeating motion indicator is obtainable by means of the performance indicator coefficient c, which is designed as c=ζ(θ(t)−θ(0)), where ζ represents the response coefficient to a joint offset; d) converting the time-varying convex quadratic programming problem, into a time-varying matrix equation by using a Lagrangian function; e) solving the time-varying matrix equation by the variable parameter convergence differential neural network; wherein the error of the time-varying matrix equation converges to zero based on the variable parameter convergence differential neural network, an error function is constructed as: ε(t)=Qy−u where ε(t) represents the error of the time-varying matrix equation, and then based on a neurodynamic method, the error is designed to converge to zero in the following way: $\frac{d\;{ɛ(t)}}{dt} = {{- \gamma}\;{\exp(t)}{\Phi\left( {ɛ(t)} \right)}}$ where γ is a parameter for adjusting a convergence rate, Φ(•) is an activation function, and the error function is substituted into the above formula to obtain a variable parameter convergence differential neural network solver, namely: Q{dot over (y)}=−{dot over (Q)}y−γ exp(t)Φ(Qy−u)+{dot over (u)} where the variable parameter convergence differential neural network solver obtains an optimal solution y* of the time-varying matrix equation, and the first n terms thereof are an optimal solution x* of the time-varying convex quadratic programming problem, i.e., an optimal solution of a joint angular velocity; f) integrating the optimal solution x* of the robotic arm at the velocity level to obtain an optimal solution θ* of a joint angle; and g) using the optimal solution θ* of the joint angle to operate the robotic arm to perform the repeating motion, so that after the end of the robotic arm completes a cycle of actions, all the joints of the robotic arm can return to an initial position.
 2. The method for performing a repeating motion of a robotic arm on the basis of a variable parameter convergence differential neural network, as claimed in claim 1, wherein the inverse kinematics equation of the robotic arm is expressed as: f(θ)=r where r is a desired track of an end of the robotic arm, and f(•) is a nonlinear equation of the joint angle of the robotic arm, and the inverse kinematics equation of the robotic arm at the velocity level is obtained by deriving the two sides of the equation with respect to time: J(θ){dot over (θ)}={dot over (r)} where J(θ)∈R^(m×n) is an m×n-dimensional matrix on the real number field, J(θ) is a Jacobian matrix of the redundant robotic arm, n is the number of the degrees of freedom of the robotic arm, and m is the number of the spatial dimensions of the track of the end of the robotic arm, and {dot over (θ)} and {dot over (r)} are respectively the derivatives of the joint angle and the track of the end of the redundant robotic arm with respect to time.
 3. The method for performing a repeating motion of a robotic arm on the basis of a variable parameter convergence differential neural network, as claimed in claim 1, wherein the formula for creating the inverse kinematics problem as a time-varying convex quadratic programming problem constrained by an equality is: $\min\left( {\frac{x^{T}{Wx}}{2} + {c^{T}x}} \right)$ s.t.  (J(θ)x = b) where x={dot over (θ)}, b={dot over (r)}, W=I represents an identity matrix, J(θ) is the Jacobian matrix of the robotic arm, and c is a performance indicator coefficient.
 4. The method for performing a repeating motion of a robotic arm on the basis of a variable parameter convergence differential neural network, as claimed in claim 1, wherein the Lagrangian function is constructed as: ${L\left( {x,\lambda,t} \right)} = {\frac{x^{T}{Wx}}{2} + {c^{T}x} + {\lambda\left( {{{J(\theta)}x} - b} \right)}}$ where λ is a Lagrangian multiplier, and the partial derivatives of the Lagrangian function respectively with respect to x and λ are obtained: $\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial x} = {{{Wx} + c + {{J^{T}(\theta)}\lambda}} = 0}$ $\frac{\partial{L\left( {x,\lambda,t} \right)}}{\partial\lambda} = {{{{J(\theta)}x} - b} = 0}$ the above system of equations can be expressed as the following time-varying matrix equation:      Qy = u ${{{where}\mspace{14mu} Q} = {\begin{bmatrix} W & {J^{T}(\theta)} \\ {J(\theta)} & 0 \end{bmatrix} \in E^{{({n + m})} \times {({n + m})}}}},{y = {\begin{bmatrix} x \\ \lambda \end{bmatrix} \in R^{({n + m})}}},{u = {\begin{bmatrix} {- c} \\ b \end{bmatrix} \in {R^{({n + m})}.}}}$ 